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Abstract 

Gravitational wave emission from the gravitational collapse of mas- 
sive stars has been studied for more than three decades. Current state 
of the art numerical investigations of collapse include those that use pro- 
genitors with more realistic angular momentum profiles, properly treat 
microphysics issues, account for general relativity, and examine non- 
axisymmetric effects in three dimensions. Such simulations predict that 
gravitational waves from various phenomena associated with gravitational 
collapse could be detectable with ground-based and space-based interfer- 
ometric observatories. 
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1 Introduction 



The field of gravitational wave (GW) astronomy will soon become a reality. The 
first generation of ground-based interferometric detectors (LIGO [1] , VIRGO [2] , 
GEO 600 [3], TAMA 300 [4]) are beginning their search for GWS. Towards 
the end of this decade, two of these detectors (LIGO, VIRGO) will undergo 
upgrades that should allow them to reach sensitivities necessary to regularly 
detect emission from astrophysical sources. A space-based interferometric de- 
tector, LISA [5], could be launched in the early part of the next decade. One 
important class of sources for these observatories is stellar gravitational col- 
lapse. This class includes the accretion induced collapse (AIC) of white dwarf 
binary components and the core collapse of massive stars (M > 8M Q ), very 
massive Population III stars (M=100-500M Q ), and supermassive stars (SMSs, 
M > 10 6 Mq). Some of these collapses result in explosions (Type II, Ib/c super- 
novae and hypernovae) and all leave behind neutron star or black hole remnants. 

Strong GWs can be emitted during a gravitational collapse/explosion and, 
following the collapse, by the resulting compact remnant [197, 142, 146, 71, 175, 
78, 79, 101]. GW emission during the collapse itself may result if the collapse 
or explosion involves aspherical bulk mass motion or convection. Rotational or 
fragmentation instabilities encountered by the collapsing star will also produces 
GWs. Neutron star remnants of collapse may emit GWs due to the growth of 
rotational or r-mode instabilities. Black hole remnants will also be sources of 
GWs if they experience accretion induced ringing. All of these phenomena have 
the potential of being detected by gravitational wave observatories because they 
involve the rapid change of dense matter distributions. 

Observation of gravitational collapse by gravitational wave detectors will 
provide unique information, complementary to that derived from electromag- 
netic and neutrino detectors. Gravitational radiation arises from the coher- 
ent superposition of mass motion. Whereas, electromagnetic emission is pro- 
duced by the incoherent superposition of radiation from electrons, atoms, and 
molecules. Thus GWs carry different kinds of information than other types of 
radiation. Furthermore, electromagnetic radiation interacts strongly with mat- 
ter and thus only gives a view of the collapse from lower density regions near 
the surface of the star and is weakened by absorption as it travels to the de- 
tector. Neutrinos can escape from much further within the collapsing star, but 
even they are scattered by the highest density regions in the core. By contrast, 
gravitational waves can propagate from the innermost parts of the stellar core 
to detectors without attenuation by intervening matter. 

The characteristics of the GW emission from gravitational collapse have been 
the subject of much study. Core collapse supcrnovae, in particular, have been 
investigated as sources of gravitational radiation for more than three decades 
(see, e.g., [166, 199, 167, 57, 151, 143, 189, 72, 141, 217, 162, 78, 79]). However, 
during this time research has produced estimates of GW strength that vary 
over orders of magnitude. This is due to the complex nature of core collapse. 
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Important theoretical and numerical issues include 

• construction of accurate progenitor models, including realistic angular mo- 
mentum distributions, 

• proper treatment of microphysics, including the use of realistic equations 
of state and neutrino transport, 

• simulation in three-dimensions to study non-axisymmetric effects, 

• inclusion of general relativistic effects, 

• inclusion of magnetic field effects, and 

• study of the effect of an envelope on core behavior. 

To date, collapse simulations generally include state of the art treatments of only 
one or two of the above physics issues (often because of numerical constraints). 
For example, those studies that include advanced microphysics have often been 
run with Newtonian gravity (and approximate evaluation of the GW emission; 
see section 2.4). A 3D, general relativistic collapse simulation, which includes 
all significant physics effects, is not feasible at present. However, good progress 
has been made on the majority of the issues listed above; the more recent work 
will be reviewed in some detail here. 

The remainder of this article is structured as follows. Each category of 
gravitational collapse will be discussed in a separate section (AIC in section 2, 
collapse of massive stars in section 3, collapse of Population III stars in sec- 
tion 4, and collapse of SMSs in section 5). Each of these sections (2, 3, 4, 
5) is divided into subsection topics: Collapse Scenario, Formation Rate, GW 
Emission Mechanisms, and Numerical Predictions of GW Emission. In the sub- 
sections on numerical predictions, the detectability of the GW emission from 
various phenomena associated with collapse is examined. In particular, the pre- 
dicted characteristics of GW emission are compared to the sensitivities of LIGO 
(for sources with frequencies of 1 to 10 4 Hz) and LISA (for sources with lower 
frequencies in the range of 10~ 4 to 1 Hz). 
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2 Accretion Induced Collapse 



2.1 Collapse Scenario 

As a white dwarf accretes matter from a binary companion, its density will 
increase. If the rate of the resulting compressional heating exceeds the white 
dwarf's cooling rate, the star's temperature will increase as well. When the 
mass of the white dwarf exceeds the Chandrasckhar stability limit, it will begin 
to collapse. 

If the temperature of the core is high enough, nuclear burning will begin and 
the stellar pressure will increase. However, electron capture will act to reduce 
both the temperature and pressure behind the burn front. The collapse out- 
come is determined by the relative strength of the nuclear burning and electron 
capture. If nuclear burning is strong enough, the white dwarf will explode as 
a Type la supernova (SN). If electron capture wins out, the collapse will con- 
tinue and result in neutron star formation. Electron capture will dominate if 
the density at which nuclear ignition occurs exceeds a critical density p cr n ■ For 
C-0 white dwarfs, p crit is in the range 6 x 10 9 - 10 10 gcm~ 3 [34]. For O-Ne- 
Mg white dwarfs, electron capture may be stronger than nuclear burning under 
most conditions (if a Rayleigh- Taylor instability does not produce a turbulent 
burn front) [157, 105]. The collapse of an O-Ne-Mg white dwarf begins when its 
central density reaches 4 x 10 9 gcm -3 . (For more details about the conditions 
under which AIC occurs, sec [157, 105, 34, 33, 129].) 

The dynamics of the collapse itself are somewhat similar to the dynamics 
of core collapse SNe (without the resulting explosive disruption of the star). 
The collapse is halted when nuclear densities are reached in the core. The core 
bounces and sends a bounce-shock outward through the star. After several 
milliseconds, the shock stalls and the remainder of the star collapses through 
the shock front (see, e.g., [211]). Less than 10 _1 Mq will likely be ejected by the 
star (due to the bounce itself or to neutrino absorption/ wind mechanisms) [75]. 

2.2 Formation Rate 

The AIC occurrence rate is difficult to determine for a number of reasons. These 
include incomplete understanding of binary star evolution and accretion pro- 
cesses [47, 116]. Another uncertainty is whether the collapse of an accreting 
Chandrasckhar mass white dwarf results in a supernovae explosion or a com- 
plete AIC (with accompanying neutron star formation). 

The AIC rate can be indirectly inferred from the observed amount of rare, 
neutron rich isotopes present in the Galaxy. These isotopes (formed via electron 
capture) are present in the portion (~ 0.1M Q ) of the outer envelope ejected by 
the star during an AIC. If all of these isotopes present in the Galaxy arc assumed 
to have originated in AICs, an upper limit of ~ 10 _5 yr _1 can be set for the 
galactic AIC rate [75] . 
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Binary population synthesis analysis can be used to determine which ac- 
creting white dwarfs will undergo AIC. The results of Yungelson and Livio [214] 
predict that the galactic AIC rate is between 8 x 1CP 7 and 8 x 10 _5 yr _1 . Thus, a 
reasonable occurrence rate can be found for an observation distance of 100 Mpc. 

2.3 GW Emission Mechanisms 

During AIC, emission of GWs will occur if the infall of matter is aspherical. GWs 
will also be produced if the collapsing star or neutron star remnant encounters 
non-axisymmetric instabilities. These include global rotational mode, r-mode, 
and fragmentation instabilities. 

Global rotational instabilities in fluids arise from non-axisymmetric modes 
e ±4m *, where m=2 is known as the "bar-mode" [193, 13]. It is convenient 
to parameterize a system's susceptibility to these modes by the stability pa- 
rameter — T rot /\W\. Here T rot is the rotational kinetic energy and W is 
the gravitational potential energy. Dynamical rotational instabilities, driven by 
Newtonian hydrodynamics and gravity, develop on the order of the rotation pe- 
riod of the object. For the uniform density, incompressible, uniformly rotating 
Maclaurin spheroids, the dynamical bar-mode instability sets in at 0d ~ 0.27. 
For differentially rotating fluids with a polytropic equation of state, numeri- 
cal simulations have determined that the stability limit 0d ~ 0.27 is valid for 
initial angular momentum distributions that are similar to those of Maclaurin 
spheroids [181, 61, 135, 103, 159, 99, 202]. If the object has an off-center den- 
sity maximum, 0d could be as low as 0.10 [201, 210, 159, 48]. Furthermore, 
recent work by Shibata, Kariono, and Eriguchi [185] suggests that 0d could be 
as low as ~ 0.01 for stars with a large degree of differential rotation. General 
relativity may enhance the dynamical bar-mode instability by slightly reducing 
0d [183, 171]. Secular rotational instabilities are driven by dissipative processes 
such as gravitational radiation reaction and viscosity. When this type of insta- 
bility arises, it develops on the timescale of the relevant dissipative mechanism, 
which can be much longer than the rotation period. The secular bar-mode 
instability limit for Maclaurin spheroids is S w 0.14. 

In rotating stars, gravitational radiation reaction drives the r-modes towards 
unstable growth [14, 74]. In hot, rapidly rotating neutron stars, this instability 
may not be suppressed by internal dissipative mechanisms (such as viscosity 
and magnetic fields) [126]. If not limited, the dimensionless amplitude a of the 
dominant (m=2) r-mode will grow to order unity within ten minutes of the 
formation of a neutron star rotating with a millisecond period. The GWs emit- 
ted carry away angular momentum, and will cause the newly formed neutron 
star to spin down over time. The spindown timescale and the strength of the 
GWs themselves are directly dependent on the maximum value a max to which 
the amplitude is allowed to grow [127, 128]. Originally, it was thought that 
otmax ~ 1- Later work indicated that a max may be > 3 [127, 190, 174, 128]. 
Some research suggests that magnetic fields, hyperon cooling, and hyperon bulk 
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viscosity may limit the growth of the r-mode instability, even in nascent neu- 
tron stars [115, 114, 164, 165, 128, 125, 87, 13] (significant uncertainties remain 
regarding the efficacy of these dissipative mechanisms). Furthermore, a recent 
study of a simple barotropic neutron star model by Arras et al. [16] suggests that 
multimode couplings could limit a max to values <C 1. If a max is indeed <C 1, 
GW emission from r-modes in collapsed remnants is likely undetectable. For 
the sake of completeness, an analysis of GW emission from r-modes (which as- 
sumes a max <~ 1) is presented in the remainder of this paper. However, because 
it is quite doubtful that a max is sizeable, r-mode sources are omitted from fig- 
ures comparing source strengths and detector sensitivities and from discussions 
of likely detectable sources in the concluding section. 

There is some numerical evidence that a collapsing star may fragment into 
two or more orbiting clumps [82]. If this does indeed occur, the orbiting frag- 
ments would be a strong GW source. 

2.4 Numerical Predictions of GW Emission 
2.4.1 Full Collapse Simulations 

The accretion-induced collapse of white dwarfs has been simulated by a num- 
ber of groups [21, 137, 211, 75]. The majority of these simulations have been 
Newtonian and have focused on mass ejection and neutrino and 7-ray emission 
during the collapse and its aftermath (note that neglecting relativistic effects 
likely introduces an error of order (v 2 /c 2 ) <~ 10% for the neutron star remnants 
of AIC; see below). The most sophisticated are those carried out by Fryer et 
al. [75], as they include realistic equations of state, neutrino transport, and 
rotating progenitors. 

As a part of their general evaluation of upper limits to GW emission from 
gravitational collapse, Fryer, Holz, and Hughes (hereafter, FHH) [78] examined 
an AIC simulation (model 3) of Fryer et al. [75]. FHH used both numerical 
and analytical techniques to estimate the peak amplitude h p k, energy Eqw, 
and frequency few of the gravitational radiation emitted during the collapse 
simulations they studied. 

For direct numerical computation of the GWs emitted in these simula- 
tions, FHH used the quadrupole approximation, valid for nearly Newtonian 
sources [140]. This approximation is standardly used to compute the GW emis- 
sion in Newtonian simulations. The reduced or traceless quadrupole moment of 
the source can be expressed as 



where i,j=l, 2, 3 are spatial indices and r—(x 2 + y 2 + z 2 ) 1 / 2 is the distance to 
the source. The two polarizations of the gravitational wave field, h + and h x , 

can be computed in terms of ly, where an overdot indicates a time derivative 




(1) 
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d/dt. The energy Eqw is a function of ly. Equation (1) can be used to 
calculate 1^- from the results of a numerical simulation by direct summation over 
the computational grid. Numerical time derivatives of 1^ can then be taken 
to compute h + , h x , and Eqw- However, successive application of numerical 
derivatives generally introduces artificial noise. Methods for computing and 

lij without taking numerical time derivatives have been developed [72, 28, 
141]. These methods recast the time derivatives of as spatial derivatives 
of hydrodynamic quantities computed in the collapse simulation (including the 
density, velocities, and gravitational potential) . Thus, instantaneous values for 

iij and 1^ can be computed on a single numerical time slice (see [141, 217] for 
details). Note that FHH define h p k as the maximum value of the rms strain 

h = {h\ + h%), where angular brackets indicate that averages have been 
taken over both wavelength and viewing angle on the sky. 

Errors resulting from the neglect of general relativistic effects (in collapse 
evolutions as a whole and in GW emission estimations like the quadrupole ap- 
proximation) are of order v 2 /c 2 . These errors are typically ~ 10% for neutron 
star remnants of AIC, < 30% for neutron star remnants of massive stellar col- 
lapse, and > 30% for black hole remnants. Neglect of general relativity in rota- 
tional collapse studies is of special concern because relativistic effects counteract 
the stabilizing effects of rotation (see section 3.4). 

Because the code used in the collapse simulations examined by FHH [75] 
was axisymmetric, their use of the numerical quadrupole approximation dis- 
cussed above does not account for GW emission that may occur due to non- 
axisymmetric mass flow. The GWs computed directly from their simulations 
come only from polar oscillations (which are significant when the mass flow 
during collapse [or explosion] is largely aspherical). 

In order to predict the GW emission produced by non-axisymmetric insta- 
bilities, FHH employed rough analytical estimates. The expressions they used 
to approximate the rms strain h and the power P = dEcw I dt of the GWs 
emitted by a star that has encountered the bar-mode instability are 

[32Gmr 2 Lu 2 

hbar = V4E-?—d- (2) 

and 

P b ar = ~mW. (3) 

45 c° 

Here m, 2r, u are the mass, length, and angular frequency of the bar and d is 
the distance to the source. FHH vary the mass m assumed to be enclosed by 
the bar (which has a corresponding length 2r) and compute the characteristics 
of the GW emission as a function of this enclosed mass. For simplicity, FHH 
assumed that a fragmentation instability will cause a star to break into two 
clumps (although more clumps could certainly be produced). Their estimates 
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for the rms strain and power radiated by the orbiting binary fragments are 



/l28G mrV 
V 5 c 4 d 



(4) 



and 




(5) 



Here m is the mass of a single fragment, 2r is the separation of the fragments, 
u) is their orbital frequency. 

For their computation of the GWs radiated via r-modes, FHH used the 
method of Ho and Lai [96] (which assumes a max =l) and calculated only the 
emission from the dominant m=2 mode. This approach is detailed in FHH. If 
the neutron star mass and initial radius are taken to be 1.4M Q and 12.53 km, 
respectively, the resulting formula for the average GW strain is 



where a is the mode amplitude and v s is the spin frequency. 

FHH's numerical quadrupole estimate of the GWs from polar oscillations in 
the AIC simulation of Fryer et al. [75] , predicts a peak dimensionless amplitude 
h pk = 5.9x 1(T 24 (for d=100 Mpc). The energy E GW = 3x 10 45 ergs is emitted at 
a frequency of few ~ 50 Hz. This amplitude is about an order of magnitude too 
small to be observed by the advanced LIGO-II detector. The sensitivity curve 
for the broadband configuration of the LIGO-II detector is shown in figure 1 
(see appendix A of [78] for details on the computation of this curve) . Note that 
the characteristic strain h is plotted along vertical axis in figure 1 (and in LISA's 
sensitivity curve, shown in figure 12). For burst sources, h = h p k- For sources 
which persist for N cycles, h — ^f~Nh p k- 

The simulation of Fryer et al. [75] considered by FHH does produce an object 
with an off-center density maximum. However, because the maximum (3 reached 
in the AIC simulation of Fryer et al. was < 0.06 (and because the degree of 
differential rotation present in the remnant was low [76]), the collapsing object 
is not likely to encounter dynamical rotational instabilities. The choice of initial 
angular momentum J for the progenitor in an AIC simulation can affect this 
outcome. Fryer et al.'s choice of J=10 49 gcm 2 s _1 is such that the period of the 
white dwarf progenitor is 10s less than the shortest observed period 30s of a 
cataclysmic variable white dwarf [117]. If higher values of J exist in accreting 
white dwarfs, bar-mode instabilities may be more likely to occur (see discussion 
of work by Liu and Lindblom below) . 

According to FHH, the remnant of this AIC simulation will be susceptible 
to r-mode growth. Assuming a max ~ 1 (which is likely not physical; see sec- 
tion 2.3), they predict Eqw > 10 52 ergs. FHH compute h(fow) for coherent 
observation of the neutron star as it spins down over the course of a year. For 



h(t) = 1.8 x 10~ 24 a 



llkHz/ \ d ) 



I v s X / 20Mpc \ 



(6) 
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Figure 1: A comparison between the GW amplitude h(f) for various sources and 
the LIGO-II sensitivity curve. See the text for details regarding the computations 
of h. The AIC sources are assumed to be located at a distance o/lOOMpc; the 
SNe sources at 10 Mpc; and the Population III sources at a luminosity distance 
of ~ 50Gpc. Secular bar-mode sources are identified with an (s), dynamical 
bar-modes with a (d). 



a neutron star located at a distance of lOOMpc, this track is always below the 
LIGO-II noise curve. The point on this track with the maximum h, which cor- 
responds to the beginning of r-mode evolution, is shown in figure 1. The track 
moves down and to the left (i.e., h and / decrease) in this figure as the r-mode 
evolution continues. 

2.4.2 Analyses of Equilibrium Representations of Collapse 

In addition to full hydrodynamics collapse simulations, many studies of gravi- 
tational collapse have used hydrostatic equilibrium models to represent stars at 
various stages in the collapse process. Some investigators use sequences of equi- 
librium models to represent snapshots of the phases of collapse (e.g., [23, 155, 
129]). Others use individual equilibrium models as initial conditions for hydro- 
dynamical simulations(e.g., [188, 159, 154, 48]). Such simulations represent the 
approximate evolution of a model beginning at some intermediate phase during 
collapse or the evolution of a collapsed remnant. These studies do not typically 
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follow the intricate details of the collapse itself. Instead, their goals include de- 
termining the stability of models against the development of non-axisymmetric 
modes and estimation of the characteristics of any resulting GW emission. 

Liu and Lindblom [129, 130] have applied this equilibrium approach to AIC. 
Their investigation began with a study of equilibrium models built to repre- 
sent neutron stars formed from AIC [129]. These neutron star models were 
created via a two step process, using a Newtonian version of Hachisu's self- 
consistent field method [83]. Hachisu's method ensures that the forces due to 
the centrifugal and gravitational potentials and the pressure are in balance in 
the equilibrium configuration. 

Their process of building the nascent neutron stars began with the construc- 
tion of rapidly rotating, pre-collapse white dwarf models. Their Models I and II 
are C-0 white dwarfs with central densities pc^lO 10 and 6 x 10 9 gcm~ 3 , respec- 
tively (recall this is the range of densities for which AIC is likely for C-0 white 
dwarfs) . Their Model III is an O-Ne-Mg white dwarf that has p c = 4 x 10 9 g cm~ 3 
(recall this is the density at which collapse is induced by electron capture). All 
three models are uniformly rotating, with the maximum allowed angular veloci- 
ties. The models' values of total angular momentum are roughly 3-4 times that 
of Fryer et al.'s [75] AIC progenitor Model 3. The realistic equation of state 
used to construct the white dwarfs is a Coulomb corrected, zero temperature, 
degenerate gas equation of state [172, 51]. 

In the second step of their process, Liu and Lindblom [129] built equilibrium 
models of the collapsed neutron stars themselves. The mass, total angular 
momentum, and specific angular momentum distribution of each neutron star 
remnant is identical to that of its white dwarf progenitor (see section 3 of [129] 
for justification of the specific angular momentum conservation assumption). 
These models were built with two different realistic neutron star equations of 
state. 

Their cold neutron star remnants had values of the stability parameter (3 
ranging from 0.23-0.26. It is interesting to compare these results with those 
of Zwerger and Muller [217]. Zwerger and Miiller performed axisymmetric hy- 
drodynamics simulations of stars with polytropic equations of state (P oc p T ). 
Their initial models were T=4/3 polytropes, representative of massive white 
dwarfs. All of their models started with ( o c =10 10 gcm -3 . Their model that was 
closest to being in uniform rotation (A1B3) had 22% less total angular momen- 
tum than Liu and Lindblom's Model I. The collapse simulations of Zwerger and 
Muller that started with model A1B3 all resulted in remnants with values of 
j3 < 0.07. Comparison of the results of these two studies could indicate that 
the equation of state may play a significant role in determining the structure of 
collapsed remnants. Or it could suggest that the assumptions employed in the 
simplified investigation of Liu and Lindblom are not fully appropriate. Zwerger 
and Muller's work will be discussed in much more detail in section 3, as it was 
performed in the context of core collapse supernovae. 

In a continuation of the work of Liu and Lindblom, Liu [130] used linearized 
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hydrodynamics to perform a stability analysis of the cold neutron star AIC 
remnants of Liu and Lindblom [129]. He found that only the remnant of the 
O-Ne-Mg white dwarf (Liu and Lindblom's Model III) developed the dynamical 
bar-mode (m=2) instability. This model had an initial /3=0.26. Note that the 
m=l mode, observed by others to be the dominate mode in unstable models 
with values of j3 much lower than 0.27 [201, 210, 159, 48], did not grow in his 
simulation. Because Liu and Lindblom's Models I and II had lower values of 
(3, Liu identified the onset of instability for neutron stars formed via AIC as 
(3 d « 0.25. 

Liu estimated the peak amplitude of the GWs emitted by the Model III 
remnant to be h p k ~ 1.4 x 10 -24 and the LIGO-II signal to noise ratio (for a 
persistent signal like that seen in the work of [154] and [36]) to be S/N < 3 
(for few ~ 450 Hz. These values are for a source located at lOOMpc. He also 
predicted that the timescale for gravitational radiation to carry away enough 
angular momentum to eliminate the bar-mode is tqw ~ 7s (~ 3 x 10 3 cycles). 
Thus, h ~ 8 x 10~ 23 . (Note this value for h is merely an upper limit as it assumes 
that the amplitude and frequency of the GWs do not change over the 7 s during 
which they are emitted. They will of course change as angular momentum is 
carried away from the object via GW emission.) Such a signal may be marginally 
detectable with LIGO-II (see figure 1). Details of the approximations on which 
these estimates are based can be found in [130]. 

Liu cautions that his results hold if the magnetic field of the proto-neutron 
star is B < 10 12 G. If the magnetic field is larger, then it may have time to 
suppress some of the neutron star's differential rotation before it cools. This 
would make bar formation less likely. Such a large field could only result if 
the white dwarf progenitor's B field was > 10 s G. Observation based estimates 
suggest that about 25% of white dwarfs in interacting close binaries (cataclysmic 
variables) are magnetic and that the field strengths for these stars are ~ 10 7 - 
3 x 10 8 G [206]. 

2.5 Going Further 

The AIC scenario is generally discussed in terms of the collapse of an accret- 
ing white dwarf to a neutron star. However, Shibata, Baumgarte, and Shapiro 
have examined the collapse of a rotating, supramassive neutron star to a black 
hole [184]. Such supramassive neutron stars (with masses greater than the 
maximum mass for a nonrotating neutron star) could be formed and pushed to 
collapse via accretion from a binary companion. They performed 3D, fully gen- 
eral relativistic hydrodynamics simulations of uniformly rotating neutron stars. 
Dynamical non-axisymmetric instabilities (such as the bar-mode) did not have 
time to grow in their simulations prior to black hole formation. Differentially ro- 
tating neutron star progenitors could have higher values of (3 than the uniformly 
rotating models used in this study and may be susceptible to non-axisymmetric 
instabilities on a shorter timescale. 
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3 Collapse of Massive Stars 



3.1 Collapse Scenario 

Stars with mass greater than ~ 8M will undergo core collapse at the end of 
their thermonuclear burning life cycles. The collapse will trigger a supernova 
explosion if the star's mass is < 40 — 5OM (for reviews, see [15, 26, 40]). Core 
collapse SNe include Types II and Ib/Ic. SNe Ib/Ic are distinguished from Type 
II SNe by the absence of hydrogen in their spectra. SNe Ib/Ic are thought to 
result from the collapse of the cores of massive stars that have lost their hydrogen 
envelopes (and possibly part of their helium envelopes) by stellar winds or by 
mass transfer. The SNe Ib/Ic progenitors that lose their outer envelopes via 
stellar winds are known as Wolf-Rayet stars and have initial masses ~ 3OM ; 
those that undergo mass loss via mass transfer in binaries have progenitor masses 
in the range 12M -18M [146]. 

It is believed that the remnant of a core collapse SN is a neutron star if the 
mass of the progenitor is less than ~ 20 — 25M Q [80, 78]. If the progenitor's 
mass is in the range 20 — 25 ~ M ~ 40 — 5OM , the entire star is not ejected 
in the SNe explosion. More than 2M will fall back onto the nascent neutron 
star and lead to black hole formation. If the progenitor's mass exceeds roughly 
40 — 5OM , then no explosion will occur and the star will collapse directly to 
a black hole. These objects are known as collapsars. However, it is unclear if 
high metallicity stars with M ~ 40 — 5OM actually end their lives in collapse 
or are prevented from doing so by stellar wind driven mass loss. [78] Note that 
the limits on the progenitor masses quoted in this paragraph (especially the 
40 — 5OM lower limit for direct black hole formation) are uncertain because the 
progenitor mass dependence of the neutrino explosion mechanism (see below) is 
unknown [88, 150]. 

The massive iron cores of SNe II/Ib/Ic progenitors are supported by both 
thermal and electron degeneracy pressures. The density and temperature of 
such a core will eventually rise, due to the build up of matter consumed by 
thermonuclear burning, to the point where electron capture and photodissocia- 
tion of nuclei begin. Dissociation lowers the photon and electron temperatures 
and thereby reduces the core's thermal support [69]. Electron capture reduces 
the electron degeneracy pressure. One or both of these processes will trigger 
the collapse of the core. The relative importance of dissociation and electron 
capture in instigating collapse is determined by the mass of the star [69]. The 
more massive the core, the bigger the role played by dissociation. 

Approximately 70% of the inner portion of the core collapses homologously 
and subsonically. The outer core collapses at supersonic speeds [69, 141]. The 
maximum velocity of the outer regions of the core reaches ~ 7 x 10 4 kms -1 . It 
takes just 1 s for an earth-sized core to collapse to a radius of 50 km [15]. 

The inward collapse of the core is halted by nuclear forces when its central 
density p c is 2-10 times the density of nuclear material [20, 19]. The core over- 
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shoots its equilibrium position and bounces. A shock wave is formed when the 
supersonically infalling outer layers hit the rebounding inner core. If the inner 
core pushes the shock outwards with energy E > 10 51 ergs (supplied by the bind- 
ing energy of the nascent neutron star), then the remainder of the star can be 
ejected in about 20 ms [15]. This so-called "prompt explosion" mechanism will 
only succeed if a very soft supra-nuclear equation of state is used in conjunc- 
tion with a relatively small core (M ~ 1.35M , derived from a very low mass 
progenitor) and a large portion of the collapse proceeds homologously [26, 146]. 
Inclusion of general relativistic effects in collapse simulations can increase the 
success of the prompt mechanism in some cases [20, 191]. 

Both dissociation of nuclei and electron capture can reduce the ejection en- 
ergy, causing the prompt mechanism to fail. The shock will then stall at a radius 
in the range 100-200km. Colgate and White [54] were the first to suggest that 
energy from neutrinos emitted by the collapsed core could revive the stalled 
shock. (See Burrows and Thompson [44] for a recent review of core collapse 
neutrino processes.) However, their simulations did not include the physics 
necessary to accurately model this "delayed explosion" mechanism. Wilson 
and collaborators were the first to perform collapse simulations with successful 
delayed ejections [32, 31, 207, 27, 209, 208]. However, their simulations and 
those of others had difficulty producing energies high enough to match observa- 
tions [52, 38, 110]. 

Observations of SN 1987a show that significant mixing occurred during this 
supernovae (see Arnett et al. [15] for a review). Such mixing can be attributed 
to nonradial motion resulting from fluid instabilities. Convective instabilities 
play a significant role in the current picture of the delayed explosion mecha- 
nism. The outer regions of the nascent neutron star are convectively unstable 
after the shock stalls (for an interval of 10- 100ms after bounce) due to the pres- 
ence of negative lepton and energy gradients [146]. This has been confirmed by 
both 2D and 3D simulations [41, 42, 110, 144, 145, 111, 146, 138, 106, 108, 161]. 
Convective motion is more effective at transporting neutrinos out of the proto- 
neutron star than is diffusion. Less than 10% of the neutrinos emitted by the 
neutron star need to be absorbed and converted to kinetic energy for the shock 
to be revived [146]. The "hot bubble" region above the surface of the neutron 
star has also been shown to be convectively unstable [52, 26, 53, 81]. Janka and 
Miiller have demonstrated that convection in this region only aids the explosion 
if the neutrino luminosity is in a narrow region [112]. Some simulations that 
include advanced neutrino transport methods have cast doubt on the ability 
of convection to ensure the success of the delayed explosion mechanism [139]. 
However, Rampp and Janka [108, 161] have recently performed collapse simu- 
lations with sophisticated neutrino physics (including Boltzmann transport and 
state of the art neutrino-matter interactions) that are very close to producing 
explosions. Their results indicate that successful explosions could be very sen- 
sitive to the details of the neutrino physics. Multi-dimensional simulations that 
include the full effects of general relativity, rotation, and an improved treatment 
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of neutrino interactions will likely be necessary to properly model the delayed 
explosion mechanism [109, 108]. 

In addition to the mixing seen in SN 1987a, observations of polarization 
in the spectra of several core collapse SNe, jets in the Cas A remnant, and 
kicks in neutron stars suggest that supernovae are inherently asphcrical (see [15, 
11, 97] and references therein). (Note that these asphcricitics could originate 
in the central explosion mechanism itself and/or the mcchanism[s] for energy 
transfer between the core and ejecta [108] .) These observations partly motivated 
the multidimensional studies of convection in the delayed explosion mechanism 
discussed above and have also spurred investigations of the role jets may play in 
SNe explosions [133, 205, 212]. Recent 3D radiation hydrodynamics simulations 
performed by Hdflich et al. [97] indicate that low velocity jets stalled inside SN 
envelopes can account for the observed asymmetries. Possible sources of such 
asymmetries include magnetic fields [22, 11], anisotropic neutrino emission, and 
convection. 

3.2 Formation Rate 

Type II/Ib/Ic supernovae are observed to occur in only spiral and irregular 
galaxies. The most thorough computation of SNe rates is that of Cappcllaro et 
al. [45]. Their sample includes 137 SNe from five different SN searches. They 
determined that the core collapse SNe rate in the Galaxy is 0.6 x 10~ 3 -1.6 x 
10~ 2 yr _1 . Thus, a reasonable occurrence rate can be found for an observation 
distance of lOMpc. (Note that a recent infrared survey estimates that the 
rate may actually be an order of magnitude higher [134].) Approximately 5 — 
40% of these SNe will leave behind black hole remnants [80]. The formation 
rate of collapsars, the massive cores that collapse to black holes without an 
accompanying SNe explosion, is unknown. This is because of the uncertainty in 
stellar wind-driven mass-loss rates [80] . 

3.3 GW Emission Mechanisms 

Gravitational radiation will be emitted during the collapse/explosion of a core 
collapse SNe due to the star's changing quadrupole moment. A rough descrip- 
tion of the possible evolution of the quadrupole moment is given in the remainder 
of this paragraph. During the first 100-250 ms of the collapse, as the core con- 
tracts and flattens, the magnitude of the quadrupole moment 1^ will increase. 
The contraction speeds up over the next 20 ms and the density distribution be- 
comes a centrally condensed torus [141]. In this phase the core's shrinking size 
dominates its increasing deformation and the magnitude of decreases. As 
the core bounces, Ijj changes rapidly due to the deceleration and rebound. If 
the bounce occurs because of nuclear pressure, its timescale will be < 1ms. If 
centrifugal forces play a role in halting the collapse, the bounce can last up to 
several ms [141]. The magnitude of lij will increase due to the core's expansion 
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after bounce. As the resulting shock moves outwards, the unshocked portion 
of the core will undergo oscillations, causing Ijj to oscillate as well. The shape 
of the core, the depth of the bounce, the bounce timescale, and the rotational 
energy of the core all strongly affect the GW emission. For further details 
see [69, 141]. 

Convectively driven inhomogeneities in the density distribution of the outer 
regions of the nascent neutron star and anisotropic neutrino emission are other 
sources of GW emission during the collapse/explosion [43, 149]. 

As discussed in the case of AIC, global rotational instabilities (such as the 
m=2 bar-mode) may develop during the collapse itself or in a neutron star 
remnant. A neutron star remnant will likely also be susceptible to the radiation 
reaction driven r-modcs. Both of these types of instabilities will emit GWs, as 
will a fragmentation instability if one occurs. See section 2.3 for further details 
regarding these instabilities. 

If the collapsed remnant is a black hole, GWs will be radiated as the infall of 
the remaining stellar matter distorts the black hole's geometry. This "ringdown 
phase" will end when gravitational radiation has dissipated all of the black hole's 
accretion-induced distortion. Zanotti, Rczzolla, and Font [215] have recently 
suggested that the torus of matter surrounding the black hole may be an even 
stronger source of GWs than the collapse itself (see also [118, 204]). When 
perturbed, such "toroidal neutron stars" may undergo regular oscillations and 
thus emit copious GWs. 

3.4 Numerical Predictions of GW Emission 
3.4.1 Historical Investigations 

The collapse of the progenitors of core collapse supernovae has been investigated 
as a source of gravitational radiation for more than three decades. In an early 
study published in 1971, Ruffini and Wheeler [166] identified mechanisms related 
to core collapse that could produce GWs and provided order of magnitude 
estimates of the characteristics of such emission. 

Quantitative computations of GW emission during the infall phase of col- 
lapse were performed by Thuan and Ostrikcr [199] and Epstein and Wag- 
oner [64, 63], who simulated the collapse of oblate dust spheroids. Thuan 
and Ostriker used Newtonian gravity and computed the emitted radiation in 
the quadrupole approximation. Epstein and Wagoner discovered that post- 
Newtonian effects prolonged the collapse and thus lowered the GW luminosity. 
Subsequently, Novikov [158] and Shapiro and Saenz [178, 167] included inter- 
nal pressure in their collapse simulations and were thus able to examine the 
GWs emitted as collapsing cores bounced at nuclear densities. The quadrupole 
GWs from the ringdown of the collapse remnant were initially investigated by 
the perturbation study of Turner and Wagoner [203] and later by Saenz and 
Shapiro [168, 169]. 
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Miiller [143] calculated the quadrupole GW emission from 2D axisymmctric 
collapse based on the Newtonian simulations of Miiller and Hillebrandt [147] 
(these simulations used a realistic equation of state and included differential 
rotation). He found that differential rotation enhanced the efficiency of the GW 
emission. 

The first fully general rclativistic investigation of stellar core collapse were 
Nakamura's 2D simulations of rotating collapse [151, 152]. However, because 
of the limits of his numerical formalism and computational resources, he was 
unable to compute the emitted gravitational radiation (the energy of this emis- 
sion is quite small compared to the rest mass energy and thus was difficult to 
extract numerically). The results of this work indicate that collapse does not 
lead to black hole formation if the parameter q = J/M 2 exceeds unity (here J 
and M are the angular momentum and gravitational mass of the remnant). 

Stark and Piran [189, 160] were the first to compute the GW emission from 
fully rclativistic collapse simulations, using the ground-breaking formalism of 
Bardeen and Piran [18]. They followed the (pressure-cut induced) collapse of 
rotating polytropes in 2D. Their work focused in part on the conditions for 
black hole formation and the nature of the resulting ringdown waveform (which 
they found could be described by the quasi-normal modes of a rotating black 
hole). In each of their simulations, less than 1% of the gravitational mass was 
converted to GW energy. 

Seidel and collaborators also studied the effects of general relativity on the 
GW emission during collapse and bounce [176, 177]. They employed a pertur- 
bative approach, valid only in the slowly rotating regime. 

The gravitational radiation from non-axisymmetric collapse was investigated 
by Detweiler and Lindblom who used a sequence of non-axisymmetric ellipsoids 
to represent the collapse evolution [57]. They found that the radiation from 
their analysis of non-axisymmetric collapse was emitted over a more narrow 
range of frequency than in previous studies of axisymmetric collapse. 

For further discussion of the first two decades of study of the GW emission 
from stellar collapse see [70]. In the remainder of section 3.4, more recent 
investigations will be discussed. 

3.4.2 Axisymmetric Simulations 

The core collapse simulations of Monchmeyer et al. began with better iron core 
models and a more realistic microphysics treatment (including a realistic equa- 
tion of state, electron capture, and a simple neutrino transport scheme) than any 
previous study of GW emission from axisymmetric stellar core collapse [141]. 
(The shortcomings of their investigation included initial models that were not 
in rotational equilibrium, an equation of state that was somewhat stiff in the 
subnuclear regime, and use of Newtonian gravity.) Each of their four models 
had a different initial angular momentum profile. The rotational energies of the 
models ranged from 0.1-.45 of the maximum possible rotational energy. 
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The collapses of three of the four models of Monchmeyer ct al. were halted by 
centrifugal forces at subnuclear densities. This type of low p c bounce had been 
predicted by Shapiro and Lightman [179] and Tohline [200] (in the context of 
the "fizzler" scenario for failed supernovae; see also [89, 90, 104]) and had been 
observed in earlier collapse simulations [148, 192]. Monchmeyer and collabo- 
rators found that a bounce caused by centrifugal forces would last for several 
ms, whereas a bounce at nuclear densities would occur in < 1 ms. They also 
determined that a subnuclear bounce produced larger amplitude oscillations in 
density and radius, with larger oscillation periods, than a bounce initiated by 
nuclear forces alone. They pointed out that these differences in timescale and 
oscillatory behavior should affect the GW signal. Therefore, the GW emission 
could indicate whether the bounce was a result of centrifugal or nuclear forces. 

Monchmeyer ct al. identified two different types of waveforms in their mod- 
els (computed using the numerical quadrupole approximation discussed in sec- 
tion 2.4). The waveforms they categorized as Type I (similar to those observed 
in previous collapse simulations [143, 72]) are distinguished by a large amplitude 
peak at bounce and subsequent damped ringdown oscillations. They noted that 
Type I signals were produced by cores that bounced at nuclear densities (or 
bounces at subnuclear densities if the cores had small ratios of radial kinetic 
to rotational kinetic energies). The quadrupole gravitational wave amplitude 
Afo for a Type I waveform is shown in figure 2 (see [195, 217] for expressions 
relating Af 2 to h) . The waveforms identified as Type II exhibit several maxima, 
which result from multiple bounces. See figure 3 for an example of a Type II 
waveform. Note that the waveforms displayed in figures 2, 3 are from the study 
of Zwerger and Mullcr [217], discussed below. 

The model of Monchmeyer ct al. that bounced due to nuclear forces had the 
highest GW amplitude of all of their models, h p k ~ 10~ 23 for a source distance 
d = lOMpc, and largest emitted energy Eqw ~ 10 47 erg. The accompanying 
power spectrum peaked in the frequency range 5 x 10 2 -10 3 Hz. 

The most extensive Newtonian survey of the parameter space of axisymmct- 
ric, rotational core collapse is that of Zwerger and Muller [217]. They simulated 
the collapse of 78 initial models with varying amounts of rotational kinetic 
energy (reflected in the initial value of the stability parameter /3j), differen- 
tial rotation, and equation of state stiffness. In order to make this large survey 
tractable, they used a simplified equation of state and did not explicitly account 
for electron capture or neutrino transport. Their initial models were constructed 
in rotational equilibrium via the method of Eriguchi and Muller [65] . The mod- 
els had a polytropic equation of state, with initial adiabatic index I\ = 4/3. 
Collapse was induced by reducing the adiabatic index to a value T r in the range 
1.28-1.325. The equation of state used during the collapse evolution had both 
polytropic and thermal contributions. 

The major result of Zwerger and Miiller's investigation was that the signal 
type of the emitted gravitational waveform in their runs was determined by the 
stiffness of the equation of state of the collapsing core (i.e., the value of T r ). In 
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Figure 2: Type I waveform (quadrupole amplitude Afo as a function of time) 
from one of Zwerger and Muller's [217] simulations of a collapsing polytrope. 
The vertical dotted line marks the time at which the first bounce occurred. (Fig- 
ure 5d of [217]; used with permission.) 

their simulations, Type I signals (as labelled by Monchmeyer et al. [141]) were 
produced by models with relatively soft equations of state T r ~ 1.31. Type II 
signals were produced by the models with stiffer equations of state (T r ~ 1.32). 
They found a smooth transition between these signal types if T r was increased 
while all other parameters were held fixed. They also observed another class 
of signal, Type III, for their models with the lowest T r =1.28. Type III wave- 
forms have a large positive peak just prior to bounce, a smaller negative peak 
just after bounce, and smaller subsequent oscillations with very short periods 
(see figure 4). Type III signals were not observed in the evolution of strongly 
differentially rotating IY=1.28 models. Their waveforms were computed with 
the same technique used in [141]. 

In contrast to the results of Monchmeyer et al. [141], in Zwerger and Muller's 
investigation, the value of p c at bounce did not determine the signal type. In- 
stead, the only effect on the waveform due to p c was a decrease in h p k in models 
that bounced at subnuclear densities. The effect of the initial value of on h p k 
was non-monotonic. For models with fa ~ 0.1, h p k increased with increasing 
0i. This is because the deformation of the core is larger for faster rotators. 
However, for models with larger (3i, h p k decreases as [3i increases. These mod- 
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Figure 3: Type II waveform (quadrupole amplitude Afo as a function of time) 
from one of Zwerger and Miiller's [217] simulations of a collapsing polytrope. 
The vertical dotted line marks the time at which bounce occurred. (Figure 5a of 
[21 7]; used with permission.) 

els bounce at subnuclear densities. Thus the resulting acceleration at bounce, 
and GW amplitude, are smaller. Zwerger and Muller found that the maximum 
value of h p k for a given sequence was reached when p c at bounce was just less 
than Pnuc- The degree of differential rotation did not have a large effect on the 
emitted waveforms computed by Zwerger and Muller. However, they did find 
that models with soft equations of state emitted stronger signals as the degree 
of differential rotation increased. 

The models of Zwerger and Muller that produced the largest GW signals 
fell into two categories: those with stiff equations of state and pi < .01; and 
those with soft equations of state, /?, > .018, and large degrees of differential 
rotation. The GW amplitudes emitted during their simulations were in the 
range 4 x 10~ 25 ~ h ~ 4 x 10~ 23 , for d=10Mpc (the model with the highest h 
is identified in figure 1). The corresponding energies ranged from 10 44 ~ E gw ~ 
10 47 erg. The peaks of their power spectra were between 500 Hz and 1 kHz. Such 
signals would fall just outside of the range of LIGO-II. 

Yamada and Sato [213] used techniques very similar to those of Zwerger and 
Muller [217] in their core collapse study. Their investigation revealed that the 
hpk for Type I signals saturated when the dimensionless angular momentum of 
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Figure 4: Type III waveform (quadrupole amplitude Afo as a function of time) 
from one of Zwerger and Miiller's [217] simulations of a collapsing polytrope. 
The vertical dotted line marks the time at which bounce occurred. (Figure 5e of 
[217]; used with permission). 

the collapsing core, q — J/(2GM/c), reached ~ 0.5. They also found that h p k 
was sensitive to the stiffness of the equation of state for densities just below 
Pnuc- The characteristics of the GW emission from their models were similar to 
those of Zwerger and Miiller. 

3.4.3 Non Axisymmetric Simulations 

The GW emission from non-axisymmetric hydrodynamics simulations of stellar 
collapse was first studied by Bonazzola and Marck [136, 30]. They used a New- 
tonian, pseudo-spectral hydrodynamics code to follow the collapse of polytropic 
models. Their simulations covered only the pre-bounce phase of the collapse. 
They found that the magnitudes of h p k in their 3D simulations were within a 
factor of two of those from equivalent 2D simulations and that the gravitational 
radiation efficiency did not depend on the equation of state. 

The first 3D hydrodynamics collapse simulations to study the GW emission 
well beyond the core bounce phase were performed by Rampp, Miiller, and 
Ruffcrt [162]. These authors started their Newtonian simulations with the only 
model (A4B5G5) of Zwerger and Miiller [217] that had a post-bounce value of 
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the stability parameter /3=0.35 that significantly exceeded 0.27 (recall this is 
the value at which the dynamical bar instability sets in for Maclaurin spheroid- 
like models). This model had the softest EOS (IY=1.28), highest $=0.04, 
and largest degree of differential rotation of all of Zwerger and Midler's models. 
The model's initial density distribution had an off-center density maximum (and 
therefore a torus-like structure). Rampp, Mullcr, and Ruffcrt evolved this model 
with a 2D hydrodynamics code until its [3 reached w 0.1. At that point, 2.5 ms 
prior to bounce, the configuration was mapped onto a 3D nested cubical grid 
structure and evolved with a 3D hydrodynamics code. 

Before the 3D simulations started, non-axisymmetric density perturbations 
were imposed to seed the growth of any non-axisymmetric modes to which the 
configuration was unstable. When the perturbation imposed was random (5% 
in magnitude), the dominant mode that arose was m=4. The growth of this 
particular mode was instigated by the cubical nature of the computational grid. 
When an m=3 perturbation was imposed (10% in magnitude), three clumps 
developed during the post-bounce evolution and produced three spiral arms. 
These arms carried mass and angular momentum away from the center of the 
core. The arms eventually merged into a bar-like structure (evidence of the 
presence of the m=2 mode). Significant non-axisymmetric structure was visible 
only within the inner 40 km of the core. Their simulations were carried out to 
<~ 14 ms after bounce. 

The amplitudes of the emitted gravitational radiation (computed in the 
quadrupole approximation) were only ^2% different from those observed in 
the 2D simulation of Zwerger and Miiller. Because of low angular resolution 
in the 3D runs, the energy emitted was only 65% of that emitted in the corre- 
sponding 2D simulation. 

The findings of Centrella et al. [48] indicate it is possible that some of the 
post-bounce configurations of Zwerger and Miiller, which have lower values of [3 
than the model studied by Rampp, Miiller, and Ruffert [162], may also be sus- 
ceptible to non-axisymmetric instabilities. Centrella et al. have performed 3D 
hydrodynamics simulations of r=1.3 polytropes to test the stability of config- 
urations with off-center density maxima (as are present in many of the models 
of Zwerger and Mullcr [217]). The simulations carried out by Centrella and 
collaborators were not full collapse simulations, but rather began with differ- 
entially rotating, equilibrium models. These simulations tracked the growth 
of any unstable non-axisymmetric modes that arose from initial 1% random 
density perturbations imposed. Their results indicate that such models can 
become dynamically unstable at values of (3 ~ 0.14. The instability observed 
had a dominant m=l mode. Centrella et al. estimate that if a stellar core of 
mass M <~ 1.4M Q and radius R ~ 200 km encountered this instability, the val- 
ues of h p k from their models would be <~ 2 x 10~ 24 -2 x 10~ 23 , for d^lOMpc. 
The frequency at which h p k occurred in their simulations was ~ 200 Hz. This 
instability would have to persist for at least <~ 15 cycles to be detected with 
LIGO-II. 
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Brown [37] carried out an investigation of the growth of non-axisymmetric 
modes in post-bounce cores that was similar in many respects to that of Rampp, 
Muller, and Ruffcrt [162]. He performed 3D hydrodynamics simulations of the 
post-bounce configurations resulting from 2D simulations of core collapse. His 
pre-collapse initial models are T = 4/3 polytropes in rotational equilibrium. 
The differential rotation laws used to construct Brown's initial models were 
motivated by the stellar evolution study of Heger, Langer, and Woosley [93]. 
The angular velocity profiles of their pre-collapse progenitors were broad and 
Gaussian-like. Brown's initial models had peak angular velocities ranging from 
0.8-2.4 times those of [93]. The model evolved by Rampp, Muller, and Ruf- 
fert [162] had much stronger differential rotation than any of Brown's models. 
To induce collapse, Brown reduced the adiabatic index of his models to r=1.28, 
the same value used by [162]. 

Brown found that (3 increased by a factor ~ 2 during his 2D collapse simu- 
lations. This is much less than the factor of ~ 9 observed in the model studied 
by Rampp, Muller, and Ruffcrt [162]. This is likely a result of the larger degree 
of differential rotation in the model of Rampp et al. 

Brown performed 3D simulations of the two most rapidly rotating of his 
post-bounce models (models £124 and £120, both of which had (3 > 0.27 after 
bounce) and of the model of Rampp et al. (which, although it starts out with 
/3=0.35, has a sustained (3 < 0.2). Brown refers to the Rampp et al. model as 
model RMR. Because Brown's models do not have off-center density maxima, 
they are not expected to be unstable to the m=l mode observed by Centrella et 
al. [48]. He imposed random 1% density perturbations at the start of all three of 
these 3D simulations (note this perturbation was of a much smaller amplitude 
than those imposed by [162]). 

Brown's simulations determined that both his most rapidly rotating model 
£124 (with post-bounce (3 > 0.35) and model RMR are unstable to growth of 
the m=2 bar-mode. However, his model £120 (with post-bounce (3 > 0.3) was 
stable. Brown observed no dominant to=3 or 4 modes growing in model RMR 
at the times at which they were seen in the simulations of Rampp et al. This 
suggests that the mode growth in their simulations was a result of the large 
perturbations they imposed. The m=2 mode begins to grow in model RMR at 
about the same time as Rampp et al. stopped their evolutions. No substantial 
m=l growth was observed. 

The results of Brown's study indicate that the overall (3 of the post-bounce 
core may not be a good diagnostic for the onset of instability. He found, as did 
Rampp, Muller, and Ruffert [162], that only the innermost portion of the core 
(with p > 10 10 gcm~ 3 ) is susceptible to the bar-mode. This is evident in the 
stability of his model £120. This model had an overall (3 > 0.3 but an inner core 
with (3i C = 0.15. Brown also observed that the (3 of the inner core does not have 
to exceed 0.27 for the model to encounter the bar-mode. Models £124 and RMR 
had [3i C ~ 0.19. He speculates that the inner cores of these later two models 
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may be bar unstable because interaction with their outer envelopes feeds the 
instability or because (3d < 0.27 for such configurations. 

Fryer and Warren [81] have recently performed the first 3D collapse simu- 
lations to follow the entire collapse through explosion. They used a smoothed 
particle hydrodynamics code, a realistic equation of state, the flux-limited dif- 
fusion approximation for neutrino transport, and Newtonian spherical gravity. 
Their initial model was nonrotating. Thus, no bar-mode instabilities could de- 
velop during their simulations. The only GW emitting mechanism present in 
their models was convection in the core. The maximum amplitude h of this 
emission, computed in the quadrupole approximation, was ~ 3 x 1CP 26 , for 
<i=10Mpc [79]. Preliminary results from rotating simulations performed with 
this same 3D code do not show the development of bar or fragmentation insta- 
bilities [79]. 

The GW emission from nonradial quasinormal mode oscillations in proto- 
neutron stars has recently been examined by Ferrari, Miniutti, and Pons [68]. 
They found that the frequencies of emission few during the first second af- 
ter formation (600-1100 Hz for the first fundamental and gravity modes) are 
significantly lower than the corresponding frequencies for cold neutron stars 
and thus reside in the bandwidths of terrestrial interferometers. However, for 
first generation interferometers to detect the GW emission from an oscillating 
proto-neutron star located at lOMpc, with a signal to noise ratio of 5, Eaw 
must be~ 10 _3 -10 _2 M©c 2 . It is unlikely that this much energy is stored in 
these modes (the collapse itself may only emit ~ 10~ 7 M Q c 2 in gravitational 
waves [59]). 

3.4.4 General Relativistic Simulations 

General relativistic effects oppose the stabilizing influence of rotation in pre- 
collapse cores. Thus stars that might be prevented from collapsing due to ro- 
tational support in the Newtonian limit, may collapse when general relativis- 
tic effects are considered. Furthermore, general relativity will cause rotating 
stars undergoing collapse to bounce at higher densities than in the Newtonian 
case [193, 217, 162, 39]. 

The full collapse simulations of Fryer and Heger [77] are the most sophis- 
ticated axisymmetric simulations from which the resultant GW emission has 
been studied [78, 79]. Fryer and Heger include the effects of general relativity, 
but assume (for the purposes of their gravity treatment only) that the mass dis- 
tribution is spherical. The GW emission from these simulations was evaluated 
with cither the quadrupole approximation or simpler estimates (see below). 

The work of Fryer and Heger [77] is an improvement over past collapse 
investigations because it starts with rotating progenitors evolved to collapse 
with a stellar evolution code (which incorporates angular momentum transport 
via an approximate diffusion scheme) [91], incorporates realistic equations of 
state and neutrino transport, and follows the collapse to late times. The values 
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of total angular momentum of the inner cores of Fryer and Hcger (0.95-1.9 x 
10 49 gcm 2 s~ 1 ) are lower than has often been assumed in studies of the GW 
emission from core collapse. Note that the total specific angular momentum of 
these core models may be lower by about a factor of 10 if magnetic fields were 
included in the evolution of the progenitors [11, 92, 94]. 

FHH's [78] numerical quadrupole estimate of the GWs from polar oscillations 
in the collapse simulations of Fryer and Heger [77] predicts a peak dimcnsionless 
amplitude h pk = 4.1 x 10~ 23 (for d=10Mpc), emitted at f GW w 20 Hz. The 
radiated energy Eqw ~ 2 x 10 44 ergs. This signal would be just out of the 
detectability range of the LIGO-II detector. 

The cores in the simulations of Fryer and Heger [77] are not compact enough 
(or rotating rapidly enough) to develop bar instabilities during the collapse and 
initial bounce phases. However, the explosion phase ejects a good deal of low 
angular momentum material along the poles in their evolutions. Therefore, 
about 1 s after the collapse, (3 becomes high enough in their models to exceed 
the secular bar instability limit. The (3 of their model with the least angu- 
lar momentum actually exceeds the dynamical bar instability limit as well (it 
contracts to a smaller radius and thus has a higher spin rate then the model 
with higher angular momentum). FHH (and [79]) compute an upper limit (via 
equation [2]) to the emitted amplitude from their dynamically unstable model of 
h ~ 3 x 10~ 22 (if coherent emission from a bar located at 10 Mpc persists for 100 
cycles). The corresponding frequency and maximum power are few ~ 10 3 Hz 
and Paw — 10 53 ergs _1 . LIGO-II should be able to detect such a signal (see 
figure 1, where FHH's upper limit to h for this dynamical bar-mode is identi- 
fied). 

As mentioned above, the proto-neutron stars of Fryer and Heger are likely 
unstable to the development of secular bar instabilities. The GW emission 
from proto-neutron stars that are secularly unstable to the bar-mode has been 
examined by Lai and Shapiro [122, 120]. Because the timescale for secular 
evolution is so long, 3d hydrodynamics simulations of the nonlinear development 
of a secular bar can be impractical. To bypass this difficulty, Lai [120] considers 
only incompressible fluids, for which there are exact solutions for (Dedekind and 
Jacobi-like) bar development. He predicts that such a bar located at lOMpc 
would emit GWs with a peak characteristic amplitude h ~ 10 , if the bar 
persists for 10 2 -10 4 cycles. The maximum few of the emitted radiation is in 
the range 10 2 — 10 3 Hz. This type of signal should be easily detected by LIGO- 

1 (although detection may require a technique like the fast chirp transform 
method of Jenet and Prince [113] due to the complicated phase evolution of the 
emission) . 

FHH predict that a fragmentation instability is unlikely to develop during 
core-collapse SNe because the cores have central density maxima (see also [79]). 
However, they do give estimates (calculated via equations [4, 5]) for the am- 
plitude, power, and frequency of the emission from such an instability: h p k ~ 

2 x 10~ 22 , P GW = 10 54 ergs-\ f GW ~ 2 x 10 3 Hz. Again, this signal would 
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fall just beyond the upper limit of LIGO-IFs frequency range. Drawing on 
analogy with fragmentation observed in star formation simulations, Davies et 
al. [56] have recently suggested that the fragmentation of collapsing massive 
stellar cores and the subsequent gravitational radiation driven coalescence of 
the fragments, may power gamma-ray bursts (see also [118]). They say that 
due to kick velocity constraints, this would be possible only for supernovae that 
leave behind black hole remnants with M ~ 12M Q (formed via the merger of the 
fragments). Under the assumption that the collapsing core splits into two frag- 
ments (each with M = O.7M ), they estimate that LIGO-II should detect the 
merger of <~ 400 core collapse fragments per year. Note that the signal strength 
would be even greater if the total mass of the black hole remnant were indeed 
12M0 (however the formation rate for these more massive objects is likely lower 
than that of lower mass cores). 

The GW emission from r-mode unstable neutron star remnants of core col- 
lapse SNe would be easily detectable if a max ~ 1 (which is likely not physical; 
see section 2.3). Multiple GW bursts will occur as material falls back onto the 
neutron star and results in repeat episodes of r-mode growth (note that a single 
r-mode episode can have multiple amplitude peaks [127]). FHH calculate that 
the characteristic amplitude of the GW emission from this r-mode evolution 
tracks from 6-1 x 10~ 22 , over a frequency range of 10 3 -10 2 Hz (see section [2.4] 
for details). They estimate the emitted energy to exceed 10 52 ergs. 

If the collapse remnant is a black hole, GWs will be emitted as the accreting 
black hole "rings down." Even with very optimistic accretion scenarios, FHH 
conclude that such radiation will be of very low amplitude and beyond the upper 
frequency reach of LIGO-II (see [78] for details). 

General relativity has been more fully accounted for in the core collapse stud- 
ies of Dimmclmeier, Font, and Muller [58, 60, 59], which build on the Newtonian, 
axisymmetric collapse simulations of Zwerger and Muller [217]. In all, they have 
followed the collapse evolution of 26 different models, with both Newtonian and 
general relativistic simulations. As in the work of Zwerger and Muller, the dif- 
ferent models are characterized by varying degrees of differential rotation, initial 
rotation rates, and adiabatic indices. They use the conformally flat metric to 
approximate the space time geometry [55] in their relativistic hydrodynamics 
simulations. This approximation gives the exact solution to Einstein's equa- 
tions in the case of spherical symmetry Thus, as long as the collapse is not 
significantly aspherical, the approximation is relatively accurate. However, the 
conformally flat condition does eliminate GW emission from the spacetime. 
Because of this, Dimmclmeier, Font, and Muller used the quadrupole approxi- 
mation to compute the characteristics of the emitted GW signal (sec [217] for 
details). 

The general relativistic simulations of Dimmclmeier ct al. showed the three 
different types of collapse evolution (and corresponding gravitational radia- 
tion signal) seen in the Newtonian simulations of Zwerger and Muller (regu- 
lar collapse-Type I signal; multiple bounce collapse-Type II signal; and rapid 
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collapse-Type III signal). However, relativistic effects sometimes led to a differ- 
ent collapse type than in the Newtonian case. This is because general relativity 
did indeed counteract the stabilizing effects of rotation and led to much higher 
bounce densities (up to 700% higher) . They found that multiple bounce collapse 
is much rarer in general relativistic simulations (occurring in only two of their 
models). When multiple bounce does occur, relativistic effects shorten the time 
interval between bounces by up to a factor of four. Movies of the simulations 
of four models from Dimmelmeier et al. [59] are shown in figures 5-8. The four 
evolutions shown include a regular collapse ( 5), a rapid collapse ( 6), a multi- 
ple bounce collapse ( 7), and a very rapidly and differentially rotating collapse 
( 8). The left frames of each movie contain the 2d evolution of the logarithmic 
density. The upper and lower right frames display the evolutions of the gravi- 
tational wave amplitude and the maximum density, respectively. These movies 
can also be viewed at [6]. 
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Figure 5: Movies from the evolution of the regular collapse model, A3B2G4, 
of Dimmelmeier et al. [59]. The left frame contains the 2d evolution of the 
logarithmic density. The upper and lower right frames display the evolutions of 
the gravitational wave amplitude and the maximum density, respectively. 
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Figure 6: Same as figure 5, but for rapid collapse model A3B2G5 of Dim- 
melmeier et al. [59]. 




Figure 7: Same as figure 5, but for multiple collapse model A2B4G1 of Dim- 
melmeier et al. [59]. 
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Figure 8: Same as figure 5, but for rapid, differentially rotating collapse model 
A4-B5G5 of Dimmelmeier et al. [59]. 

Dimmelmeier et al. found that models for which the collapse type was the 
same in both Newtonian and relativistic simulations, had lower GW amplitudes 
h p k in the relativistic case. This is because the Newtonian models were less 
compact at bounce and thus had material with higher densities and velocities 
at larger radii. Both higher and lower values of h p k were observed in models for 
which the collapse type changed. Overall, the range of h p k (4 x 10~ 24 -3 x 1CT 23 , 
for a source located at 10 Mpc) seen in the relativistic simulations was quite close 
to the corresponding Newtonian range. The average Eqw was somewhat higher 
in the relativistic case (1.5 x 10 47 ergs compared to the Newtonian value of 6.4 x 
10 46 ergs). The overall range of GW frequencies observed in their relativistic 
simulations (60-1000 Hz) was close to the Newtonian range. They did note that 
relativistic effects always caused the characteristic frequency of emission, few, 
to increase (up to five-fold) . For most of their models, this increase in few was 
not accompanied by an increase in h p k ■ This means that relativistic effects could 
decrease the detectability of GW signals from some core collapses. However, the 
GW emission from the models of Dimmelmeier et al. could be detected by the 
first generation of ground-based interferometric detectors if the sources were 
fortuitously located in the Local Group of galaxies. A catalog containing the 
signals and spectra of the GW emission from all of their models can be found 
at [7]. 

Fully general relativistic collapse simulations (i.e., without the conformally 
flat approximation) have been performed by Shibata [182]. He used an axisym- 
metric code that solves the Einstein equations in Cartesian coordinates and the 
hydrodynamics equations in cylindrical coordinates. The use of the Cartesian 
grid eliminates the presence of singularities and allows for stable, long-duration 
axisymmetric simulations [12]. The focus of this work was the effect of rota- 
tion on the criteria for prompt black hole formation. Shibata found that if the 
parameter q = J/M 2 is less than 0.5, black hole formation occurred for rest 
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masses slightly greater than the maximum mass of spherical stars. However, 
for 0.5 < q < 1, the maximum stable rest mass is increased by ~ 70-80%. The 
results are only weakly dependent on the initial rotation profile. Shibata did 
not compute the GW emission in his collapse simulations. 

The new general relativistic hydrodynamics simulations of Zanotti, Rezzolla, 
and Font [215] suggest that a torus of neutron star matter surrounding a black 
hole remnant may be a stronger source of GWs than the collapse itself. They 
used a high resolution shock-capturing hydrodynamics method in conjunction 
with a static (Schwarzschild) spacetime to follow the evolution of "toroidal neu- 
tron stars." Their results indicate that if a toroidal neutron star (with constant 
specific angular momentum) is perturbed, it could undergo regular oscillations. 
They estimate that the resulting GW emission would have a characteristic am- 
plitude h c ranging from 6 x 10~ 24 -5 x 10~ 23 , for ratios of torus mass to black hole 
mass in the range 0.1-0.5. (These amplitude values arc likely underestimated 
because the simulations of Zanotti et al. are axisymmetric.) The corresponding 
frequency of emission is few ~ 200 Hz. The values of h c and few quoted here 
are for a source located at lOMpc. This emission would be just outside the 
range of LIGO-II (sec figure 1). Further numerical investigations, which study 
tori with non-constant angular momenta and include the effects of self-gravity 
and black hole rotation, are needed to confirm these predictions. Movies from 
the simulations of Zanotti et al. can be viewed at [8]. 

Magnetized tori around rapidly spinning black holes (formed via either core 
collapse or neutron star-black hole coalescence) have recently been examined in 
the theoretical study of van Putten and Levinson [204] . They find that such a 
torus-black hole system can exist in a suspended state of accretion if the ratio 
of poloidal magnetic field energy to kinetic energy EsjE^ is less than 0.1. They 
estimate that ~ 10% of the spin energy of the black hole will be converted 
to gravitational radiation energy through multipole mass moment instabilities 
that develop in the torus. If a magnetized torus-black hole system located at 
lOMpc is observed for 2 x 10 4 rotation periods, the characteristic amplitude of 
the GW emission is ~ 6 x 10~ 20 . It is possible that this emission could take 
place at several frequencies. Observations of x-ray lines from gamma-ray bursts 
(which are possibly produced by these types of systems) could constrain these 
frequencies by providing information regarding the angular velocities of the tori 
(preliminary estimates from observations suggest few ~ 500 Hz, placing the 
radiation into a range detectable by LIGO-I). 

3.4.5 Simulations of Convective Instabilities 

Convectively driven inhomogeneities in the density distribution of the outer 
regions of the nascent neutron star and anisotropic neutrino emission are other 
sources of GW emission during the collapse/explosion [43, 146]. GW emission 
from these processes results from small-scale asphericities, unlike the large- 
scale motions responsible for GW emission from aspherical collapse and non- 
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axisymmetric global instabilities. Note that Rayleigh- Taylor instabilities also 
induce time-dependent quadrupolc moments at composition interfaces in the 
stellar envelope. However, the resultant GW emission is too weak to be detected 
because the Rayleigh- Taylor instabilities occur at very large radii [146]. 

Burrows and Hayes [43] have suggested that the large proper motions ob- 
served in pulsars may result from asymmetrical core collapse SNe. They have 
computed the resultant GW emission under the assumption that such asymme- 
tries were caused by convection and anisotropic neutrino emission during the 
collapse and explosion. Hydrodynamic simulations and theoretical investiga- 
tions suggest that asymmetries present in the star prior to collapse (in part due 
to convection during silicon and oxygen burning) will be amplified during the 
collapse [24, 121]. Based on these arguments, Burrows and Hayes induced a 
mass dipole anisotropy of < 0.1% in the pre-collapse core used in their "ex- 
ploratory" collapse simulation. The initial quadrupole anisotropy present in 
their model was 0.012. They imposed no aspherical perturbations on the initial 
velocities. The collapse was followed in ID inside a radius of 15 km and in 2D 
(with azimuthal symmetry) to larger radii. To reduce the computational burden 
of this exploratory simulation, they artificially hardened the neutrino spectrum 
to drive the explosion at a faster rate (note that this practice was observed to 
minimize the GW emission in the simulations of Miiller and Janka [149], see 
below). 

The initial asymmetries were observed to grow in Burrows and Hayes' cal- 
culation. The nascent neutron star developed a recoil speed of <~ 530 km s -1 . 
Aspherical mass motion was responsible for the vast majority of the recoil. Only 
16% of the recoil velocity could be attributed to anisotropic neutrino emission. 
The gravitational waveform from this simulation (including separate matter and 
neutrino contributions) is show in Figure 9. The contributions of matter motion 
and neutrino emission to the GW amplitude h were of opposite sign. Anisotropic 
neutrino radiation was basically the sole contributor to h for the first 20 ms, af- 
ter which the contribution from mass motion became significant. The peak 
amplitude calculated was h p k ~ 3x 10~ 24 , for a source located at lOMpc. Most 
(<~ 70%) of the emitted Eqw ~ 2.0 x 10 45 ergs came from mass motion (because 
the mass quadrupole was changing more rapidly as the bounce occurred). The 
frequency of the GWs emitted was in the range 10 < few < 100 Hz. Burrows 
and Hayes predicted LIGO-II should be able to detect this emission, with a 
signal-to-noise ratio of 10, from a core collapse that fortuitously happened to 
occur in our Galaxy (with <i=10kpc). 

The study of Nazin and Postnov [153] predicts a lower limit for Eqw emitted 
during an asymmetric core collapse SNe (where such asymmetries could be 
induced by both aspherical mass motion and neutrino emission). They assume 
that observed pulsar kicks are solely due to asymmetric collapse. They suggest 
that the energy associated with the kick (Mv 2 /2, where M and v are the mass 
and velocity of the neutron star) can be set as a lower limit for Eqw (which can 
be computed without having to know the mechanism behind the asymmetric 
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Figure 9: The gravitational waveform (including separate matter and neutrino 
contributions) from the collapse simulations of Burrows and Hayes [43] . The 
curves plot the product of the gravitational wave amplitude and distance to the 
source as a function of time. (Figure 3 of [43]; used with permission.) 

collapse). From observed pulsar proper motions, they estimate the degree of 
asymmetry e present in the collapse and the corresponding characteristic GW 
amplitude (h oc y/e). This amplitude is 3 x 10~ 25 for a source located at lOMpc 
and emitting at few = 1 kHz. 

Miiller and Janka performed both 2D and 3D simulations of convective insta- 
bilities in the proto-neutron star and hot bubble regions during the first second 
of the explosion phase of a Type II SNe [149]. They numerically computed the 
GW emission from the convection induced aspherical mass motion and neutrino 
emission in the quadrupole approximation (for details, see section 3 of their 
paper). 

For typical iron core masses, the convectively unstable region in the proto- 
neutron star extends over the inner 0.7-1.20M© of the core mass (this corre- 
sponds to a radial range of ~ 10-50 km). The convection in this region, which 
begins approximately 10-20 ms after the shock forms and may last for w 20 ms- 
1 s, is caused by unstable gradients in entropy and/or lepton number resulting 
from the stalling of the prompt shock and deleptonization outside the neutri- 
nosphcrc. Miiller and Janka's simulations of convection in this region began 
with the ID, non-rotating, 12 ms post-bounce model of Hillebrandt[95]. This 
model included general relativistic corrections that had to be relaxed away prior 
to the start of the Newtonian simulations. Neutrino transport was neglected in 
these runs (see section 2.1 of [149] for justification); however, a sophisticated 
equation of state was utilized. Figure 10 shows the evolution of the temperature 
and density distributions in the 2D simulation of Miiller and Janka. 

The peak GW amplitude resulting from convective mass motions in these 
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Figure 10: Convective instabilities inside the proto-neutron star in the 2D simu- 
lation of Mutter and Janka [14-9]- The evolutions of the temperature (left panels) 
and logarithmic density ( right panels ) distributions are shown for the radial re- 
gion 15-95 km. The upper and lower panels correspond to times 12 and 21ms, 
respectively, after the start of the simulation. The temperature values range 
from 2.5 x 10 10 to 1.8 x 10 11 K. The values of the logarithm of the density range 
from 10.5 to 13.3 gem -3 . The temperature and density both increase as the col- 
ors change from blue to green, yellow, and red. (Figure 7 of [149]; used with 
permission.) 

simulations of the proto-neutron star was w 3 x 10~ 24 in 2D and w 2 x 10~ 25 
in 3D, for c?=10Mpc. The emitted energy was 9.8 x 10 44 ergs in 2D and 1.3 x 
10 42 ergs in 3D. The power spectrum peaked at frequencies of 200-600 Hz in 
2D and 100-200 Hz in 3D. Such signals would not be detectable with LIGO-II. 
The reasons for the differences between the 2D and 3D results include smaller 
convective elements and less under- and overshooting in 3D. The relatively low 
angular resolution of the 3D simulations may have also played a role. The 
quadrupole gravitational wave amplitude A^q from the 2D simulation is shown 
in the upper left panel of figure 11 (see [217, 195] for expressions relating Af 2 
to h). 

Convection in the hot bubble region between the shock and ncutrinospherc 
arises because of an unstable entropy gradient resulting from neutrino heating. 



32 



This unstable region extends over the inner mass range 1.25-1 AQMq (corre- 
sponding to a radial range of w 100-1000 km). Convection in the hot bubble 
begins s=s 50-80 ms after shock formation and lasts from w 100-500 ms. Only 
2D simulations were performed in this case. These runs started with a 25 ms 
post-bounce model provided by Bruenn. A simple neutrino transport scheme 
was used in the runs and an imposed neutrino flux was located inside the neu- 
trinosphere. Due to computational constraints, the computational domain did 
not include the entire convectively unstable region inside the proto-neutron star 
(thus this set of simulations only accurately models the convection in the hot 
bubble region, not in the proto-neutron star). 

The peak GW amplitude resulting from these 2D simulations of convective 
mass motions in the hot bubble region was h p k «5x 10~ 25 , for c?=10Mpc. The 
emitted energy was ~ 2 x 10 42 ergs. The energy spectrum peaked at frequencies 
of 50-200 Hz. As the explosion energy was increased (by increasing the imposed 
neutrino flux), the violent convective motions turn into simple rapid expansion. 
The resultant frequencies drop to few ~ 10 Hz. The amplitude of such a signal 
would be too low to be detectable with LIGO-II. 

Miiller and Janka estimated the GW emission from the convection induced 
anisotropic neutrino radiation in their simulations (see [149] for details). They 
could not directly compute the characteristics of the GW emission from neu- 
trinos because their proto-neutron star simulations did not include neutrino 
transport and the neutrino transport in their hot bubble simulations was only 
ID. They found that the amplitude of the GWs emitted can be a factor of 
5-10 higher than the GW amplitudes resulting from convective mass motion. 
However, their results indicate that the Eqw coming from anisotropic neutrino 
emission is only a small fraction of that emitted by convective mass motion. 

Note that the qualitative results (including h and Eqw magnitudes) of Bur- 
rows and Hayes' [43] "exploratory" 2D calculation are in rough agreement with 
the more detailed 2D simulations of Miiller and Janka . However, both investiga- 
tions have shortcomings. Necessary improvements would include more realistic 
initial models (which include rotation and are produced with multi-dimensional 
collapse simulations to better estimate initial asphericities), use of sophisticated 
neutrino transport schemes, inclusion of general relativistic effects, and higher 
resolution 3D runs. 
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Figure 1 1 : Quadrupole amplitudes Af 2 [cm] from convective instabilities in vari- 
ous models of [149]. The upper left panel is the amplitude from a 2D simulation 
of proto-neutron star convection. The other three panels are amplitudes from 
2D simulations of hot bubble convection. The imposed neutrino flux in the hot 
bubble simulations increases from the top right model through the bottom right 
model. (Figure 18 of [149]; used with permission.) 
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3.5 Going Further 



The background of GW emission from a population of core-collapse SNe at 
cosmological distances may be detectable by LIGO-II, according to Ferrari, 
Matarrese, and Schneider [67]. They determined the SNe rate as a function 
of redshift using observations to determine the evolution of the star formation 
rate. Only collapses that lead to black hole formation were considered. This 
simplified the study because the GW emission from such collapses is generally a 
function of just the black hole mass and angular momentum. They found that 
the stochastic background from these sources is not continuous and suggest 
that this could be used to optimize detection strategies. The maximum GW 
spectral strain amplitude they computed was in the range 10~ 28 -1CP 27 Hz, at 
frequencies of a few times 10 2 Hz. Such a signal may be detected by a pair of 
LIGO-II detectors. 
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4 Collapse of Population III Stars 



4.1 Collapse Scenario 

The first generation of stars to form in the early universe are known as Pop- 
ulation III stars (formed at redshifts z ~ 5). Theoretical and computational 
evidence suggests that Population III stars may have had masses ~ 1OOM [10, 
9, 82]. Because these massive stars contained no metals, it was possible for them 
to form directly and then evolve with very low stellar winds and thus very little 
mass loss. If the mass of a nonrotating Population III star is ~ 260M Q , its fate 
is to collapse directly to a black hole at the end of its life [82]. If rotational 
support prevents the star from direct collapse to a black hole, explosive ther- 
monuclear burning will cause the star to undergo a giant hypernova explosion. 
Prior to black hole formation, the rotating, collapsed core will have a mass of 
50-70M Q and a radius of 1000 2000 km. Note that because these massive stars 
evolve so quickly (in a few million years [17]), the events associated with their 
deaths will take place at roughly the redshifts of their births. 

4.2 Formation Rate 

The formation rate of Population III stars can be indirectly estimated from 
the re-ionization fraction of the early universe, which was re-ionized by light 
from these stars [98, 46]. Using estimates of the ultraviolet light emission from 
Population III stars, their ionization efficiency, and the re-ionization fraction of 
the early universe, one can determine that about 0.01%-1% of the universe's 
baryonic matter was found in these very massive stars. This corresponds to ~ 
10 4 -10 7 Population III stars in a lO n M0 galaxy and thus a collapse rate that is 
~ 10~ 3 yr _1 . Thus, a reasonable occurrence rate can be found for an observation 
(luminosity) distance of ~ 50Gpc (which corresponds to a redshift of z=5, in 
the cosmology used by [78]). Uncertainties in the assumptions used make this 
formation rate uncertain to within a few orders of magnitude. [9, 82, 10, 78] It 
is unknown how many Population III stars will collapse directly to black holes. 

4.3 GW Emission Mechanisms 

The GW emission mechanisms related to the collapse of Population III stars are 
a subset of those discussed in the sections on AIC and SNe/collapsars. These 
include asphcrical collapse, global rotational and fragmentation instabilities that 
may arise during the collapse/explosion and in the collpase remnant (prior to 
black hole formation), and the "ring-down" of the remnant black hole. 
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4.4 Numerical Predictions of GW Emission 

The GW emission from the collapse of Population III stars has recently been in- 
vestigated by Fryer and collaborators (Fryer, Woosley, and Hegcr [82], FHH [78], 
and Fryer, Holz, Hughes, and Warren [79]). The collapse simulations of Fryer, 
Woosley, and Heger again started with rotating collapse progenitors that had 
been evolved with a stellar evolution code [91]. The initial models used by 
the evolution code were in rigid rotation with a surface ratio of centrifugal to 
gravitational forces of 20% (this ratio is seen in current observations of O stars) . 

The results of Fryer, Woosley, and Hegcr suggest that the collapse remnant 
(prior to black hole formation) is susceptible to the development of a secular bar- 
mode instability. However, at z > 5, the GW emission would be redshiftcd out 
of LIGO-H's frequency range. At z—5, h p k = 8 x 10~ 23 , with a corresponding 
frequency of 10 Hz [78, 79]. Even if such a signal persists for a hundred cycles, 
it would likely be undetectable by LIGO-II. Note that these signal strengths 
are orders of magnitude lower than the qualitative estimates of signal strength 
given in Carr, Bond, and Arnett [46]. 

LIGO-II may be able to detect the GW emission from binary clumps formed 
via a fragmentation instability. If such a signal is emitted at z—5 and persists 
for 10 cycles, h would be ~ 10~ 22 , over a frequency range of 10-100 Hz [78, 79]. 
The likelihood of the development of a fragmentation instability is diminished 
by the fact that the off-center density maxima present in the simulations of 
Fryer, Woosley, and Heger are not very pronounced [76]. 

The "ring-down" of the black hole remnant will likely be strong because 
Fryer, Woosley, and Heger observe a high accretion rate after collapse. FHH 
estimate that for a source located at z=20, the GWs would be redshiftcd out 
of LIGO-H's bandwidth. However, for a source at z—5, h p k ~6x 10~ 23 and 
the frequency range is 20-70 Hz. This signal may be marginally detectable with 
LIGO-II (see figure 1). 
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5 Collapse of Supermassive Stars 



5.1 Collapse Scenario 

There is a large body of observational evidence that supermassive black holes 
(SMBHs, M ~ 10 6 Mq) exist in the centers of many, if not most galaxies (see, 
e.g., the reviews of Rees [163] and Macchetto [132]). The masses of SMBHs in 
the centers of more than 45 galaxies have been estimated from observations [66] 
and there are more than 30 galaxies in which the presence of a SMBH has been 
confirmed [119]. 

One of the possible formation mechanisms for SMBHs involves the gravi- 
tational collapse of supermassive stars (SMSs). The timescale for this forma- 
tion channel is short enough to account for the presence of SMBHs at red- 
shifts z > 6 [107]. Supermassive stars may contract directly out of the pri- 
mordial gas, if radiation and/or magnetic field pressure prevent fragmenta- 
tion [86, 62, 85, 131, 35, 9]. Alternatively, they may build up from fragments 
of stellar collisions in clusters [173, 25]. Supermassive stars are radiation domi- 
nated, isentropic and convective [181, 216, 131]. Thus, they are well represented 
by an n~3 polytrope. If the star's mass exceeds 10 6 M Q , nuclear burning and 
electron/positron annihilation are not important. 

After formation, an SMS will evolve through a phase of quasistationary 
cooling and contraction. If the SMS is rotating when it forms, conservation 
of angular momentum requires that it spins up as it contracts. There are two 
possible evolutionary regimes for a cooling SMS. The path taken by an SMS 
depends on the strength of its viscosity and magnetic fields and on the nature 
of its angular momentum distribution. 

In the first regime, viscosity or magnetic fields are strong enough to en- 
force uniform rotation throughout the star as it contracts. Baumgarte and 
Shapiro [23] have studied the evolution of a uniformly rotating SMS up to the 
onset of relativistic instability. They demonstrated that a uniformly rotating, 
cooling SMS will eventually spin up to its mass shedding limit. The mass shed- 
ding limit is encountered when matter at the star's equator rotates with the 
Keplerian velocity. The limit can be represented as fished = (T/\W\) s hed- In 
this case, fished — 9 x 10~ 3 . The star will then evolve along a mass shedding 
sequence, losing both mass and angular momentum. It will eventually contract 
to the onset of relativistic instability [102, 49, 50, 181, 107]. 

Baumgarte and Shapiro used both a second-order, post-Newtonian approx- 
imation and a fully general relativistic numerical code to determine that the 
onset of relativistic instability occurs at a ratio of R/M ~ 450, where R is the 
star's radius and G—c—1 in the remainder of this section. Note that a second- 
order, post-Newtonian approximation was needed because rotation stabilizes 
the destabilizing role of nonlinear gravity at the first post-Newtonian level. If 
the mass of the star exceeds 10 6 M Q , the star will then collapse and possibly 
form a SMBH. If the star is less massive, nuclear reactions may lead to explosion 
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instead of collapse. 

The major result of Baumgarte and Shapiro's work is that the universal 
values of the following ratios exist for the critical configuration at the onset of 
relativistic instability: T/|W|, R/M, and J/M 2 . These ratios are completely 
independent of the mass of the star or its prior evolution. Because uniformly 
rotating SMSs will begin to collapse from a universal configuration, the subse- 
quent collapse and the resulting gravitational waveform will be unique. 

In the opposite evolutionary regime, neither viscosity nor magnetic fields 
are strong enough to enforce uniform rotation throughout the cooling SMS as 
it contracts. In this case it has been shown that the angular momentum distri- 
bution is conserved on cylinders during contraction [29] . Because viscosity and 
magnetic fields are weak, there is no means of redistributing angular momentum 
in the star. So even if the star starts out rotating uniformly, it cannot remain 
so. 

The star will then rotate differentially as it cools and contracts. In this 
case, the subsequent evolution depends on the star's initial angular momentum 
distribution, which is largely unknown. One possible outcome is that the star 
will spin up to mass shedding (at a different value of [3 s hed than a uniformly 
rotating star) and then follow an evolutionary path that may be similar to that 
described by Baumgarte and Shapiro [23]. The alternative outcome is that 
the star will encounter the dynamical bar instability prior to reaching the mass 
shedding limit. New and Shapiro [155, 156] have demonstrated that a bar-mode 
phase is likely to be encountered by differentially rotating SMSs with a wide 
range of initial angular momentum distributions. This mode will transport mass 
and angular momentum outward and thus may hasten the onset of collapse. 

5.2 Formation Rate 

An estimate of the rate of the collapse of SMSs can be derived from the quasar 
luminosity function. Hachnelt [84] has used the quasar luminosity function to 
compute the rate of GW bursts from supermassive black holes, assuming that 
each quasar emits one such burst during its lifetime (and that each quasar is 
a supermassive black hole). If it is assumed that each of these bursts is due 
to the formation of a supermassive black hole via the collapse of a SMS, then 
Hachnelt's rate estimates can be used as estimates of the rate of SMS collapse. 
This rate is likely an overestimate of the SMS collapse rate because many SMBHs 
may have been formed via merger. Haehnelt predicts that the integrated event 
rate through redshift z=4.5 ranges from ~ 10~ 6 yr -1 for M=10 8 M Q objects to 
~ lyr -1 for M=10 6 M Q objects. Thus, as in the case of Population III stars, a 
reasonable occurrence rate can be found for an observation (luminosity) distance 
of 50 Gpc. 
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5.3 GW Emission Mechanisms 



The GW emission mechanisms related to the collapse of SMSs are a subset 
of those discussed in the sections on AIC, SNe/collapsars, and Population III 
stellar collapse. These include aspherical collapse, global rotational and frag- 
mentation instabilities that may arise during the collapse/explosion and in the 
collapsed remnant (prior to black hole formation) , and the "ring-down" of the 
remnant black hole. 



5.4 Numerical Predictions of GW Emission 

The outcome of SMS collapse can only be determined with numerical, relativistic 
3D hydrodynamics simulations. 

Until recently, such simulations had only been published for nearly spherical 
collapse. The spherical simulations of Shapiro and Teukolsky [180] produced 
collapse evolutions that were nearly homologous. In this case, the collapse time 
Tcoii is roughly the free-fall time at the horizon 

Tcou = (J^) ' = 14s[M/lO 6 M ]-\ (7) 

The peak GW frequency f GW =r~ oU is then 10~ 2 Hz, if the mass of the star is 
10 6 M Q . This is in the middle of LISA's frequency band of 10~ 4 -1 Hz [198, 73]. 

The amplitude h of this burst signal can be roughly estimated in terms of 
the star's quadrupole moment 

2M 2 
h e- 



Rd 

< elxlO~ ls [M/lO 6 M ][d/5OGpc}- 1 . (8) 

Here d is the distance to the star and e ~ ?Y|H^| is a measure of the star's 
deviation from spherical symmetry. In this case, e will be much less than one 
near the horizon, since the collapse is nearly spherical. 

There are two possible aspherical collapse outcomes that can be discussed. 
The first outcome is direct collapse to a SMBH. In this case, e will be on the 
order of one near the horizon. Thus, according to equation 8, the peak amplitude 
of the GW burst signal will be h pk ~ 1 x 10- 18 [M/10 6 M Q ][d/50 Gpc]" 1 . 

Alternatively, the star may encounter the dynamical bar mode instability 
prior to complete collapse. Baumgarte and Shapiro [23] have estimated that a 
uniformly rotating SMS will reach (3 ~ 0.27 when R/M=15. The frequency of 
the quasiperiodic gravitational radiation emitted by the bar can be estimated 
in terms of its rotation frequency to be 



fgw — ^fbar 



( 



GM 



1/2 



V R 3 

= 2 x 1O- 3 Hz[M/1O 6 M ]-\ (9) 
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when R/M=15. The corresponding h p k, again estimated in terms of the star's 
quadrupole moment, is 

2M 2 

hpk < 



Rd 

< 1 x HT 19 [M/lO 6 M ][d/5OGpc] _1 . (10) 

The LISA sensitivity curve is shown in figure 12 (see [100] for details on the 
computation of this curve; a mission time of 3yrs has been assumed). The GW 
signal from this dynamical bar-mode could be detected with LISA. 




0.0001 0.0010 0.0100 0.1000 1.0000 

f (Hz) 

Figure 12: A comparison between the GW amplitude h(f) for various sources 
and the LISA noise curve. See the text for details regarding the computations 
of h. The SMS sources are assumed to be located at a luminosity distance of 
50Gpc. The bar-mode source is a dynamical bar-mode. 

Shibata and Shapiro [186] have recently published a fully general relativis- 
tic, axisymmetric simulation of the collapse of a rapidly, rigidly rotating SMS. 
They found that the collapse remained homologous during the early part of the 
evolution. An apparent horizon does appear in their simulation, indicating the 
formation of a black hole. Because of the symmetry condition used in their run, 
non-axisymmetric instabilities were unable to develop. 

The collapse of a uniformly rotating SMS has been investigated with post- 
Newtonian hydrodynamics, in 3+1 dimensions, by Saijo, Baumgarte, Shapiro, 
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and Shibata [170]. Their numerical scheme used a post-Newtonian approxima- 
tion to the Einstein equations, but solved the fully relativistic hydrodynamics 
equations. Their initial model was an n—S polytrope. 

The results of Saijo et al. indicate that the collapse of a uniformly rotat- 
ing SMS is coherent (i.e., no fragmentation instability develops). The collapse 
evolution of density contours from their model is shown in figure 13. Although 
the work of Baumgarte and Shapiro [23] suggests that a bar instability should 
develop prior to BH formation, no bar development was observed by Saijo et 
al. They use the quadrupole approximation to estimate a mean GW amplitude 
from the collapse itself: h = 4 x 10~ 21 , for a 1O 6 M0 star located at a distance 
of 50 Gpc. Their estimate for few at the time of BH formation is 3 x 10 -3 Hz. 
This signal would be detectable with LISA (see figure 12). 




Figure 13: Meridional plane density contours from the SMS collapse simulation 
of Saijo, Baumgarte, Shapiro, and Shibata [170]. The contour lines denote 
densities p=p c x d^^ l ^ le \ where p c is the central density. The frames are 
plotted at (t/t D , Pc , d)=(a)(5.0628 x 10~ 4 , 8.254 x 10~ 9 , 10~ 7 J, (b) (2.50259, 
1.225 x 10~ 4 , 10~ 5 j; fcj(2.05360, 8.328 x 10" 3 , 5.585 x lO" 7 ), (d)(2. 50405, 
3.425 x 10 -2 , 1.357 x 10 J, respectively. Here t, tr>, and Mq are the time, 
dynamical time ( —yj B? e jM , where R e is the initial equatorial radius and M is 
the mass), and rest mass. (Figure 15 of [170]; used with permission.) 

Saijo et al. also consider the GW emission from the ringdown of the black 
hole remnant. For the l=m=2 quasi-normal mode of a Kerr black hole with 
a/M=0.9, they estimate the characteristic frequency and amplitude of the emis- 
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sion to be few ~ 2 x 1(T 2 Hz and h ~ 1 x 10 _20 [(Ai?Gn//M)/10 _4 ] 1 / 2 , for an 
M = 10 6 M Q source located at a luminosity distance of 50Gpc (see [123, 196, 
187] for details). Here AEqw/M is the radiated energy efficiency and may be 
~ 7 x 10~ 4 [189]. This GW signal is within LISA's range of sensitivity (see 
figure 12). 
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6 Summary 



It is hoped that as gravitational collapse simulations become more sophisticated, 
the historically widely varying estimates of the magnitude of GW emission from 
collapse may start to converge. Steady progress in this field has been made in 
the last decade. Some researchers have begun to use progenitor models produced 
with stellar evolution codes, which thus have more realistic angular momentum 
profiles, as starting points for collapse simulations [78]. This reduces the need for 
collapse studies that include large surveys of the angular momentum parameter 
space. Other progress made in the numerical study of collapse includes the use 
of realistic equations of state [78] , advanced neutrino transport and interaction 
schemes [108, 124, 161, 194], and the performance of 3D Newtonian [162, 37, 81] 
and improved 2D general relativistic simulations [182]. 

There is still much work to be done toward the goal of self-consistent, 3D 
general relativistic collapse simulations. Accurate progenitor modelling and col- 
lapse simulations must include the effects of magnetic fields, as they can signifi- 
cantly alter the amount of angular momentum and differential rotation present 
in collapsing stars. Many of the more advanced studies, which include proper 
microphysics treatment and/or general relativistic effects, have been limited to 
axisymmetry. Full 3D simulations are necessary to compute the characteristics 
of the GW emission from non-axisymmetric collapse phenomena. Furthermore, 
simulations that follow both the collapse and the evolution of the collapsed 
remnant are necessary to consistently predict GW emission. One benefit of 
long duration simulations is that they will facilitate the investigation of the ef- 
fects of the envelope on any instabilities that develop in the collapsing core or 
remnant. Of course, lengthy 3D simulations are computationally intensive. This 
burden may be reduced by the use of advanced numerical techniques, including 
adaptive mesh refinement and parallel algorithms. 

The current numerical simulations of gravitational collapse indicate that 
interferometric observatories could detect GWs emitted by some collapse phe- 
nomena. LIGO-I may be able to detect GWs from secular bar-mode instabilities 
in core-collapse SNe [120] and magnetized tori surrounding black hole collapse 
remnants [204]. LIGO-II could observe GWs from dynamical bar-mode insta- 
bilities in AIC [130] and core-collapse SNe [78], and possibly from the fragmen- 
tation of very massive SNe cores that merge to form BHs [56] . LISA should be 
able to detect the collapse (and any bar-mode instabilities that develop during 
the collapse) of SMSs [23] and the ringdown of black hole remnants of collapsed 
Population III stars [78] and SMSs [170]. These observations will provide unique 
information about gravitational collapse and its associated progenitors and rem- 
nants. 
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